Conjugated donor-acceptor chromophores in solution: non-linearity at work 
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I. INTRODUCTION 

Conjugated donor-acceptor (DA) chromophores, also 
called push-pull chromophores, are molecules constituted 
by an electron-donor (D) and an electron-acceptor (A) 
group connected by a 7r-conjugated backbone. These 
molecules are widely investigated in several, apparently 
unrelated fields: they are common solvation probescl, 
are molecules of choice -for second order non-linear op- 
tics (NLO) applicatkinsa, and are useful model systems 
for electron transfeiO. All these applications exploit the 
presence of a low-lying excited state, characterized by a 
different electronic distribution from the ground state. 

Reactions involving charge migrations are ubiquitous 
in all fields of chemistry and biochemistry, so that under- 
standing their mechanism is of special interesta. These 
reactions often occur in solution, and are strongly influ- 
enced by the medium itself, suggesting a strong coupling 
between solvent and solute. Nowadays, advances in laser 
technology allow to study the coupled reactant and sol- 
vent dynamicsQ. In the typical experiment, a (series of) 
short pulse (s) of light is used to perturb the electronic 
charge distribution of the solute; the system is then in- 
terrogated by a retarded probe pulse. Polar conjugated 
chromophores, due to their strongly solvent-dependent 
properties, are molecules of choice for this kind of exper- 
iments, offering as an additional bonus also the oppor- 
tunity to investigate the kinetics, of electron transfer in 
the absence of reactant diffusiontl. A proper interpreta- 
tion of both steady-state and time-resolved experiments 
is the key to gain important physical information. How- 
ever, standard interpretative schemes are based on linear 
perturbative treatments of the solute-solvent interaction, 
so that they disregard one of the most characteristic fea- 
tures of push-pull chromophores, i.e. their large non- 
linearity. We have already proposed a model that, ac- 
counting for the intrinsic non-linearity of the electronic 
system, is able to rationalize steady-state electronic and 
vibrational spectra of polar chromophores in solution. In 
this paper the same model is extended to discuss time- 
resolved experiments. 



II. A MODEL FOR NON-LINEAR 
INTERACTIONS 

Since the properties of push-pull chromophores are 
dominated by the lowest excited state, a two-state model 
describes the essential physics governing the low-energy 
behavior of these molecules. The electronic structure of 
push-pull chromophores can be described as resonating 
between the fully neutral, \DA), and the zwitterionic, 
\D^A~), states, separated by an energy 2zo and mixed 
up by the hopping integral, — V2nj. The ionicity oper- 
ator, p, measures the weight of \D^A~) in the ground- 
or excited-state. The expectation value of this operator, 
p = (/)), corresponding to the polarity of the molecule 
in the relevant state, only depends on the ratio zq/I, so 
that the electronic problem is fully described by a single 
parameter. We allow for the coupling of the electronic 
system with a molecular vibration, Q, and a solvation 
coordinate, q, by assigning the two basis states two har- 
monic potential energy surfaces (PES), with equal cur- 
vature but displaced minima along the two coordinates. 
The resulting Hamiltonian, describing these slow degrees 
of freedom and their coupling to electrons isQ: 

n' = ^n^Q^ - y^nQp + ^UJ^q^ - V^^ujqp (1) 

where the first two terms account for the internal vibra- 
tion, and the last two for the solvation coordinate. f2 and 
uj are the frequencies corresponding to Q and q, respec- 
tively; Csp and for are the energies gained by \D^A^) due 
to the relaxation of Q and q, respectively, i.e. the small- 
polaron binding energy and the solvent reorganization 
energy. 

The linear dependence of the energy separation be- 
tween the two basis states on Q and q coordinates orig- 
inates a non-linear dependence of p on the two coordi- 
nates, making the exact potential energy, V{Q,q), an- 
harmonic. This anharmonicity, originating from the cou- 
pling of electrons to slow degrees of freedom, has many 
important consequences: in particular it is responsiWe 
for a large amplification of static hyperpolarizabilitiesQ. 



III. SPECTROSCOPIC PROPERTIES 
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Absorption and fluorescence spectra are vertical pro- 
cesses, occurring at different (Q,g). The dependence of 
p on slow coordinates (both Q and q) then leads to non- 
trivial effects on absorption and fluorescence frequen- 
cies and band-shapes. The vertical excited state reached 
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upon photon absorption has a different polarity with re- 
spect to the ground state. Immediately after absorption, 
the slow degrees of freedom readjust in response to the 
new charge distribution of the solute. But, as long as 
slow degrees of freedom relax, the solute molecule itself 
feels a new surrounding and in turn readjusts its polar- 
ity. The non-linear, self-consistent relaxation problem 
can be solved exactly in our simple picture, to calcu- 
late the equilibrium polarity of the fully relaxed excited 
stat^. Fluorescence is once more a vertical process that, 
starting from the relaxed excited state, leads to the or- 
thogonal ground state. Our model then predicts different 
Huang-Rhys factors for absorption and fluorescence, and 
explains the observation of narrower fluoreseence than 
absorption bands for push-pull chromophoresB. 

The solvent orientational coordinate describes a very 
slow (actually overdamped) motion and is responsible 
for inhomogeneous broadening effects in optical spectra. 
This broadening can be modeled in teems of a statistical 
distribution of q conformations, 'w{qf^. Since the chro- 
mophore readjusts its polarity in response to the solvent 
configuration, one ends up with a distribution of chro- 
mophore polarities, w{p). Broadened spectra are then 
calculated as a sum of spectra corresponding to different 
molecular polarities. 

The self-consistent nature of the interactions implies 
that slow degrees of freedom are in turn affected by the 
coupling to electrons. The most apparent effects can be 
found in the dependence of vibrational properties on the 
chromophore charge distribution and hence on the sol- 
vent polarity. In agreement with experimental data, we 
predict solvent-dependent vibrational frequencies and in- 
homogeneous broadening of infrared and Raman bands 
in polar solventsllB. Inhomogeneous broadening, affecting 
both electronic and vibrational states in polar solvents, 
rationalizes an unusual phenomenon of dispecsion of res- 
onant Raman bands with the excitation linetil, observed 
for the dye phenol blue dissolved in polar solvental3. 

The dependence of molecular properties and hence of 
absorption and emission band-shapes on the configura- 
tion of slow coordinates (in particular the solvation co- 
ordinate) is the key to understand some interesting phe- 
nomena observed in time-resolved spectra. Recent pump- 
probe and femtosecond hole-burning measurements_.c 
push-pull chromophores dissolved in polar solventsE 
have in fact revealed the appearance of so-called tempo- 
rary isosbestic points (TIPs) in transient gain and differ- 
ential absorption spectra, respectively. This means that 
spectra collected at different time delays seem to cross 
in an isosbestic point, whose-position actually depends 
on the chosen time windowlU. The physical origin of 
TIPs has not been-junderstood so far. The only temp- 
tative explanationll3 invokes the active role of two (or 
more) excited states. This explanation is not adequate, 
as proved by the fact that TIPs do appear in transient 
spectra collected for molecules dissolved in polar solvents 
only, |Wjhile this feature is not observed in non-dipolar sol- 
ventsta. The appearance of TIPs in transient spectra is 
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then related to solvation dynamics. 

In typical time-resolved experiments, the relaxation of 
internal molecular vibrations is in general completed af- 
ter the fast few hundreds of femtosecond following the ex- 
citations, so that subsequent dynamics can be ascribed 
to the relaxation of the solvation coordinate toward the 
new equilibrium configuration. Then, in the framework 
of our model, we can calculate the time evolution of 
transient spectra based on the temporal evolution of the 
probability distribution w(p), which, in turn, can be ob- 
tained from the evolution of w{q) on the relevant PES. 
In the phase space of coordinate and momentum {q,p) 
the evolution of w{q,p) can be described by the Fokker- 
Planck equation. In the overdamped regime (friction co- 
efficient, 7, greater than lu), relevant to the solvation 
coordinate, the Fokker-Planck equation reduces to-Jthe 
simpler Smoluchowski equation, in the g-space onlyEj: 



dw{q) 



w{q) 



d'^V dVdw{q) 



kT 



d'^w{q) 



(2) 



dq"^ dq dq dq^ 

where V is the anharmonic potential energy for the rele- 
vant state, k the Boltzmann constant, T the temperature. 
An example of the temporal evolution of w{q) on one of 
the PES is given in Fig. 1(a) from t = 0, when an out- 
of-equilibrium population is created, to i — > oo, when 
equilibrium is reached. In Fig. 1(b) the corresponding 
evolution of w(p) is reported: due to anharmonicity and 
to the non-linear relation between p and g, the distribu- 
tion is asymmetric and evolves in width and shape. 

Transient properties only depend on 7 and lo via the 
ratio r — ^ /uP'. In order not to add any adjustable pa- 
rameter to the original model, we fix r to the longitudinal 
relaxation time of the pure solvent, as obtained from lit- 
erature data. In Fig. 2(a), we show transient emission 
spectra, calculated for the parameters reported in the 
caption. The evolution of the band-shape, due to the de- 
pendence of electronic properties on the configuration of 
slow variables, together with the red-shift, is responsible 
for the appearance of a TIP. These data reproduce the 
gain spectra measured for the polyene PA2 dissolved in 
dioxanalj. The adopted r value coincides with the ac- 
cepted value of the longitudinal relaxation time of diox- 
anellj. Fig. 2(b) refers to a different experiment. In par- 
ticular, to simulate femtosecond hole-burning spectra of 
coumarin 102 (C102) dissolved in CHaCNIiJ, we have cal- 
culated transient differential absorption spectra from an 
out-of-equilibrium g-distribution in the ground state, as 
obtained by a pump-dump sequence. The spectra, cal- 
culated for thepHiodcl parameters that fit steady-state 
spectra of C102ll3, and for r — 0.3 ps (the accejjited value 
for the longitudinal relaxation time of CH3CNEJ) are in 
good agreement with experimental data. Specifically, the 
blue-shift of the transient absorption is accompanied by 
an evolution of the absorption band-shape, so that spec- 
tra collected in a narrow temporal window all cross in a 
unique point. This point is not truly isosbestic and, in 
agreement with experimental data, we calculate a slow 
blue-shift of its position at increasing time. 
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IV. CONCLUSIONS 

The proposed model describes the spectral properties 
of conjugated donor-acceptor chromophores in solution 
in terms of two electronic states coupled to internal vi- 
brations and to solvent degrees of freedom. The model is 
simple enough to allow for a full exploitation of the non- 
linearity of the electronic system, that manifests itself in 
a self-consistent interplay between electronic, vibrational 
and solvation degrees of freedom. Signatures of this non- 
linearity are recognized in several typical spectral fea- 
tures that have not explanation in standard approaches. 
The model is semiempirical in nature: the few model 
parameters are fixed against experimental data. Then 
the applicability of the model itself and the reliability 
of the model parameters for a specific molecule must be 
confirmed via the comparison with a large body of ex- 
perimental data, much wider than the pool of data used 
to fix the parameters. In this respect we underline that 
the scope of our model is extremely wide, ranging from 
steady-state to time-resolved spectra, from electronic to 
vibrational data. Based on available experimental data, 
the model has been validated for the phenol blue dye, 
whose steady-state electronic-and vibrational spectra are 
extensively discussed in Ref.til. A second interesting case 
study is offered by C102, whose time-resolved electronic 
spectra have been discussed in the previous Section. For 
this molecule, solvation effects have also been obser-«ed 
in steady-state and time-resolved vibrational spectr 
also these effects can be rationalized within our modelE^ . 
Linear electron-phonon and solute-solvent interactions 
originate, in materials with non-linear responses, large 
and non-trivial effects that cannot be understood based 
on perturbative treatments. 
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FIG. 1. (a) Evolution of the probability distribution w{q) over the represented PES. (b) Corresponding evolution of the 
probability distribution w{p). 
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(b) 
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FIG. 2. (a) Transient emission spectra calculated for zo = 0.5 cV, V^t — 1 cV, Q = 0.2 eV, esp = 0.22 cV, tor ~ 0.4 cV, 
T = 1.7 ps, T = 300 K. Time delays are displayed on the figure, time increasing in the direction of the arrows, (b) Transient 
differential absorption spectra calculated for zo = 1 eV, \/2t = 1.2 eV, Q, = 0.2 eV, esp = 0.3 eV, Coi = 0.35 eV, r = 0.3 ps, 
T = 300 K. Upper panel and lower panel refer to different time windows, as labeled on the figure. Time increases in the 
direction of the arrows. 
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